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We combine the shear-transformation-zone (STZ) theory of amorphous plasticity with Edwards' 
statistical theory of granular materials to describe shear flow in a disordered system of thermalized 
hard spheres. The equations of motion for this system are developed within a statistical thermo- 
dynamic framework analogous to that which has been used in the analysis of molecular glasses. 
For hard spheres, the system volume V replaces the internal energy U as a function of entropy 
S in conventional statistical mechanics. In place of the effective temperature, the compactivity 
X — dV/dS characterizes the internal state of disorder. We derive the STZ equations of motion 
for a granular material accordingly, and predict the strain rate as a function of the ratio of the 
shear stress to the pressure for different values of a dimensionless, temperature-like variable near a 
jamming transition. We use a simplified version of our theory to interpret numerical simulations 
by Haxton, Schmiedeberg and Liu, and in this way are able to obtain useful insights about internal 
rate factors and relations between jamming and glass transitions. 



I. INTRODUCTION 

In this paper, we combine the basic elements of the 
shear-transformation-zone (STZ) theory of amorphous 
molecular plasticity [l], [|| with Edwards' statistical the- 
ory of granular materials 0-0| to construct a theory of 
shear flow in a noncrystalline system of thermalized hard 
spheres 

The shear-transformation-zone (STZ) theory has been 
successful in accounting for a range of glassy, nonequilib- 
rium phenomena, such as the formation of shear bands 
Q, stress-strain curves of metallic glasses in constant 
strain-rate experiments |9j , and f req uency-dependent vis- 
coelastic response functions 11]. Most recently, one 
of us (JSL) has proposed that the STZ's play the role 
of dynamic heterogeneities in producing both Stokes- 
Einstcin violations and stretched-exponential relaxation 

The essence of the STZ theory of plasticity is the as- 
sumption that irreversible molecular rearrangements oc- 
cur at isolated flow defects, i.e. STZ's. In deforming sys- 
tems, the STZ's are not fixed structural features. Rather, 
they are activated fluctuations that appear and disappear 
in response to thermal or mechanically generated noise. 
The plastic strain rate in this theory is proportional to 
the STZ density, which is determined by a Boltzmann- 
like factor of the form exp (— ez/0 c g), where ez is the 
STZ formation energy, and e ff = dUc/dSc is the effec- 
tive temperature in energy units. Here, Uc is the config- 
urational potential energy and Sc is the configurational 
entropy 9 e g is a measure of the configurational disorder. 

Sheared, dense, granular materials exhibit features, 
such as dynamic yield stresses and jamming transitions, 
that are similar to those of amorphous molecular sys- 
tems. However, these systems differ fundamentally from 
each other because the energies of interaction between 
ideally hard grains are undefined. For a granular ma- 
terial composed of frictionless hard spheres, there is no 
potential energy of interaction between grains. Concepts 
such as the effective temperature # e ff and the formation 



energy ez cannot be carried over directly from theories of 
molecular plasticity. Nevertheless, granular materials do 
carry entropy, and granular rearrangements under shear 
must be governed by the second law of thermodynamics 
just as molecular rearrangements are in an elastoplastic 
solid. How, then, are we to formulate an STZ theory for 
granular materials composed of perfectly hard grains? 

Edwards and coworkers 0-0| have argued that the to- 
tal volume V occupied by a granular system plays the 
role of the configurational energy Uc in the statistical 
theory of granular materials. In general, like Uc, V is 
a function of the configurational entropy Sc plus other 
internal state variables to be specified below. The com- 
pactivity, 



X = 



dV 
dSc~ 



(1.1) 



is the analog of the effective temperature. Our main 
purpose here is to extend this analogy to nonequilibrium 
situations. 

This paper is structured as follows. We start in Sec. II 
by writing out the first and second laws of thermodynam- 
ics for a system of hard spheres in contact with a thermal 
reservoir and driven by external forces. Then, in Sees. Ill 
- V, we derive equations of motion for the STZ's and for 
the compactivity X. In Sec. VI, we use our theory to 
interpret the results of molecular dynamics simulations 
of bidisperse hard spheres by Haxton, Schmiedeberg and 
Liu [13, [HI, hereafter referred to as HSL. An especially 
important part of the HSL results is Haxton's use of a 
compressibility identity to measure an effective tempera- 
ture proportional to X 1A\ , and thereby provide informa- 
tion regarding internal rate factors and relations between 
jamming and glass transitions. We conclude in Sec. VII 
with remarks about the implications of these results. 
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II. FIRST AND SECOND LAWS OF 
THERMODYNAMICS 

Consider a noncrystalline system of hard grains whose 
kinetic temperature 9 is fixed by contact with a thermal 
reservoir. For simplicity, assume that this reservoir has 
a large enough heat capacity that its temperature does 
not change when energy flows between it and the gran- 
ular subsystem. Let Ut denote the total energy of this 
system, including the granular kinetic energies and the 
energy stored in the reservoir. If the grains interact only 
via contact forces, they have no configurational potential 
energy and, therefore, no such energy is included in Ut- 

Suppose that this system is driven in simple (not pure) 
shear by a shear stress s and a pressure p. The first law 
of thermodynamics for this system is 



Ut = V s^j — pV 

= Vs 7 — pXSc — p 



/ dV 



A Q , (2.1) 



3c 



where 7 is the shear rate, Sc is the granular configu- 
rational entropy introduced above in Eq. (jl.ll) . and the 
A Q are internal variables that specify the configurational 
state of the granular subsystem. 

Let St denote the entropy of the reservoir plus the 
(quantitatively negligible) entropy of the kinetic degrees 
of freedom of the grains. Then 



Ut — St, 



and 



pXSc = Vsj — p 



dV 



(2.2) 



A a - 0S T - (2.3) 



So 



The second law of thermodynamics requires that the to- 
tal entropy be a non-decreasing function of time: 



S = S c + S T > 



(2.4) 



Substituting Eq. (|2.3[) for Sc into the second law above, 
and using the fact that each individually variable term 
in the resulting inequality must be non- negative [H, [l6[ , 
we arrive at the second-law constraints 



w = vs 7 - p 



/ dV 



A Q >0; 



(2.5) 



\dA a/ Sc 

(pX - 6)S T > 0. (2.6) 



The dissipation rate W, as defined in [l5|, is the dif- 
ference between the rate at which inelastic work is done 
on the configurational subsystem and the rate at which 
energy is stored in the internal degrees of freedom. The 
second constraint implies that pX — 9 and St must carry 
the same sign if they are nonzero, so that 



6S T = -K 



1-^1 -a 



(2.7) 



where /C is a non- negative thermal transport coefficient. 
It is already clear from this analysis that pX plays the 
role of a temperature, p X approaches 9 in an equilibrat- 
ing system; and a heat flux Q flows from the granular 
subsystem into the reservoir when the two subsystems 
are not in thermodynamic equilibrium with each other. 



III. STZ EQUATIONS OF MOTION 

The following discussion is almost, but not quite, the 
same as that which has appeared in the STZ literature. Q 
We repeat it here because several of the differences are 
important. 

As in the case of molecular systems, we assume that 
the STZ's are two-state flow defects. We further as- 
sume that only a single species of STZ is statistically 
relevant at any given temperature, pressure, or packing 
fraction. In other words, for the steady-state behavior of 
interest here, we do not invoke the explicit distribution 
over STZ transition rates that was needed in theories of 
frequency-dependent viscoelasticity [l(J [H| or stretched- 
exponential relaxation [l2j]. However, some elements of 
those theories may be relevant here. 

Following 0, Hi, we suppose for simplicity that STZ's 
can be classified as "plus" and "minus" according to their 
orientations relative to the applied shear stress. We let 
N + and iV_ denote the number of STZs in each of the 
two orientations, and let 



N + +N_ iV + 
A = ; m = 



N_ 



N 



N+ + N- 



(3.1) 



denote the density and orientational bias of STZ's, where 
./V is the number of grains. 

Let vz denote the excess volume per STZ. Then the 
total volume V is 



V = NAv z + Vl(Si) 

= NA.v z + Vi(S c - S z {Km)), 



(3.2) 



where V± and Si are the volume and entropy of all config- 
urational degrees of freedom of the granular system not 
associated with STZ's, and Sz is the entropy associated 
with the STZ's. Then 111 



where 



S z (A,m) = NS (A) + NAif)(m) (3.3) 



5 (A) = -AlnA + A; (3.4) 
ip(m) — In 2 — — (1 + m) ln(l + m) 

--(l-m)ln(l-m). (3.5) 

The STZ equation of motion for iV + and iV_ is gener- 
ally written in the form: 

tN± = TZ(±s)N t -TZ(ts)N ± +T (^N eq - N ± ) . (3.6) 
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The corresponding strain rate is 

7 = |^ [H(s)N- K(~s)N + ] , 



(3.7) 



where, because we are describing simple rather than pure 
shear, we define the volume of the plastic core of an STZ 
to be 2vo- 

Each term in Eq. (|3.6[) needs some interpretation. On 
the left-hand side, r is a time scale that ordinarily has 
been chosen to be consistent with the underlying micro- 
scale dynamics, for example, a molecular vibration pe- 
riod. The analogous choice for the hard-sphere system is 
the inertial time scale r = ^Jm/pa, where to is the aver- 
age mass of a sphere, and a 3 is its average volume. This 
is the choice of time scale made by HSL in reporting their 
data; thus, we adopt it here. For relatively small packing 
fractions <fi = Na 3 /V, it can be shown that the inertial 
r is proportional to the average time between sphere- 
sphere collisions multiplied by a dimensionless function 
of <p] thus, the inertial time scale is a natural analog of 
the molecular vibration period. At larger values of (f>, 
however, where the system becomes jammed, we expect 
that most of the spheres are in close contact with each 
other, and that their rearrangement rates — just as in 
molecular glasses — are determined by collective motions 
that are much slower than r _1 . 

On the right-hand side of Eq. (|3.6[) . the first two 
terms containing lZ(±s) are the rates (in units of t -1 ) 
at which the STZ's are making forward and backward 
transitions. These rates describe much of the basic 
physics of this class of systems; they are discussed in 
more detail throughout this paper. The second two terms 
are the rates of STZ creation and annihilation. These 
are fluctuation-activated processes, expressed here in the 
form of a detailed-balance relation in which N eq is the 
steady-state, total number of STZ's. 

r/r is an attempt frequency consisting of additive ther- 
mal and mechanical parts: 

f = p + r. (3.8) 

The quantity p is best understood as a dimensionless, 
thermal noise strength. Well above the jamming transi- 
tion, where r is an accurate approximation for the ele- 
mentary time scale, p should be of the order of unity. As 
we approach the jamming transition, either by decreas- 
ing or increasing p or (f>, p decreases rapidly, becoming 
unmeasureably small below a transition point. When 
p = 0, the system is fully jammed in the sense that con- 
figurational rearrangements can occur only in response to 
sufficiently large driving forces. The rate p/r, multiplied 
by an activation factor of the form exp (— const. /0), is our 
analog of the "a" relaxation rate t" 1 . A first-principles 
calculation of r a is the central unsolved problem of glass 
physics. We do not attempt to solve it here; instead, we 
deduce values of p from the HSL data. 

In analogy to p/r, the quantity T/t is the contribution 
to the attempt frequency in Eq. (|3.6p due to mechanically 



generated noise. Here, we do have a well defined prescrip- 
tion for computing T from the rate of entropy generation, 
as shown below in Eq. (|4.ip . 

In terms of the intensive variables A and m, the STZ 
equations of motion become: 



tA 



f(A eq - A); 



(3.9) 

tto = 2 C(s)(T(s) — to) — f to — r— to; (3.10) 
r 7 = 2e AC(s)(T(s)-m), (3.11) 
where e = Nvq/V and A eq = N eq /N. We also define 



C(s) = ~(K(s)+K(-s)); 



and 



T(s) = 



n(s)-n(-a) 
n(s) + Ti(-s)' 



(3.12) 



(3.13) 



At this point in the development, the second law 
of thermodynamics provides useful constraints on the 
ingredients of the preceding equations. Substituting 
Eqs. (EU), (I3TTU)) and (j3~TT|) into Eq. for the dis- 

sipation rate, we find 



T — = — 1 p X A TO— 

N dm 
+2 A C(s) (T(s) - to) (v q s + P X 



dm 



-pT(A eq - A) 



vz + X I In A — il)(m) + to — — 

dm 



(3.14) 



The second-law constraint, W > 0, implies that each 
of the three terms in Eq. A3. 14)) must be non-negative. 
The first term automatically satisfies this requirement 
because, from Eq. (|3.5[) . we have 



di\) 
dm 



1 



In 



1 



1 — TO 



tanlT^TO) (3.15) 



so that the product —m(dip/dm) is automatically non- 
negative. 

An especially important result comes from the second 
term, for which 

(T(s) - to) [v s + PX^Q > 0. (3.16) 

The two factors on the left-hand side must be mono- 
tonically increasing functions of s that change sign at 
the same point for arbitrary values of to. According to 
Eq. (|3.15l) . this is possible only if 



T(s) = tanh 



vqs 
pX 



(3.17) 



The non-negativity constraint on the third term in 
Eq. (|3.14|) can be written in the form 



dF 
9A 



(A eq - A) > Q 



(3.18) 
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where F is a free energy given by 



-F(A, m) = p 



v z A + XS a (A) - XA [ tp(m) - m-^- 

dm 



(3.19) 

A eq must be the value of A at which dF/dA changes sign, 
so that 



A eq = exp 



v z , ,, , dip 
— + tphn) — m- — 
X am 



2 exp 



vz_ 
X 



(3.20) 

Thus, the STZ density in this non-equilibrium situation 
is given by a Boltzmann-like expression in which the com- 
pactivity plays the role of the temperature. 



IV. QUASISTATIONARY RELATIONS 

At this point in our analysis, we specialize to quasis- 
tationary or steady-state situations for which we can set 
A = m = 0, implying specifically that A = A eq . We start 
by assuming that Pechenik's hypothesis [HI, [ljj remains 
valid for a sheared granular material; that is, that the 
mechanical noise strength T is proportional to the heat 
production per STZ. In steady flow, all of the work done 
on the system is dissipated as heat; therefore the rate 
of energy dissipation per unit volume is 7s. To convert 
this rate into a noise strength with dimensions of inverse 
time, we multiply by the volume per noise source, i.e. the 
volume per STZ, v/A eq , and divide by an energy conve- 
niently written in the form eovso- Here, v = V/N, and 
so is a system-specific parameter with the dimensions of 
stress. The resulting expression for F is 



r 



S"fT 



2s 



e s A eq s 



C(s)(T(s) 



(4.1) 



With this result, the stationary version of Eq. (|3.10p 
reads 



2C(s)(T(s)-m) [ 1- — 



mp = 0, 



(4.2) 



which is satisfied by m — m e9 (s), where 

s , p 



m eq {s) 



£0 
2s 



1 



s 2C(s) 



1 



-T(») 



2C(s) 

In particular, when p = 0, we find 



4-T(s) 
«o 



T(s), if (s/s )T(s) < 1; 
sq/s, if (s/s )T(s) > 1. 



(4.3) 



(4.4) 



Thus an exchange of stability occurs in the same manner 
as it did in molecular systems, with the low-temperature 
yield stress being the solution of the equation 



SyT(sy) = Sy tanli 



vps y 
pX 



so- 



(4.5) 



If the temperature- like quantity p X is small in compar- 
ison with ^oSo) then s„ w so. 

Finally, the steady-state version of Eq. (|3.11[) for the 
strain rate becomes 



r 7 = 4e e- 1 / x C(s) 



tanh 



VQS 

vzPX 



m eq {s) 



(4.6) 

Here, we have introduced the inertial number q as a di- 
mensionless measure of the strain rate. (E.g. see [2(|.) 
We also have introduced the dimensionless compactivity 
X = X/v z . 



KINEMATIC EQUATIONS FOR THE 
COMPACTIVITY 



The dimensionless compactivity \ is t ne analog of the 
effective temperature in molecular systems and, as such, 
is a measure of the system's state of structural disorder. 
Referring to Eq. (|2.7|) . we see that the natural definition 
of a thermal temperature comparable to x is 

a 

(5.1) 



pvz 



which is almost the same as the dimensionless tempera- 
ture §h = 0/pa 3, defined by HSL. 

Below the jamming transition, i.e. at low tempera- 
tures or high pressures where p = 0, and where long- 
lived structural rearrangements can be induced only by 
externally driven deformation, we expect there to be a 
direct relation between the steady-state shear rate and 
the compactivity, say, x = x(?)- When the shear rate 
is much smaller than any relevant relaxation rate in the 
system, then dimensional analysis requires that x(q) be a 
q-independent constant, say xo- However, when the shear 
rate becomes comparable to internal relaxation rates, i.e., 
when the system is being "stirred" rapidly on its intrinsic 
time scales, its disorder increases, and x(q) becomes an 
increasing function of q. The simplest way to describe 
this situation is to write ^asa function of the ratio q/qo , 
where qo/r is an internal rate relevant to shear relax- 
ation. One obvious possibility is that qo is proportional 
to 72.(0), which is a measure of the rate at which STZ's 
spontaneously undergo shear transitions between their 
two states. This rate becomes slow in jammed systems, 
but it generally remains nonzero even when p = 0. 

In recent publications 0, [22|, it has been assumed 
that the large-g relation between q/qo and x nas a Ar- 
rhenius form, q/qo ~ exp(— A/jQ for x 3> Xo- Both \ 
and A are volumes measured in units of vz ■ If vz is the 
only volume in the system relevant to configurational re- 
arrangements, then we expect that A is of the order of 
unity. A relation of this form was discovered in numerical 
simulations by Haxton and Liu [l8j . 

To interpolate between the sm all-q and large-g behav- 
iors of x(<?); it was proposed in [l9j that x(q) be writ- 
ten phenomenologically in the form of a Vogel-Fulcher- 
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Tamann (VFT) expression for a "viscosity" q 1 as a func- 
tion of the temperature x : 



1 1 

- = — exp 
1 <7o 



A 

X 



with, 



a e ff(x) 



Xi 



X - Xo 



exp 



X - Xo 
XA ^ Xo 



(5.2) 



(5.3) 



Thus, x — > xo m the limit of small strain rates, and 
X — > oo as q — > go- The exponential cutoff in Eq. (|5.3p 
is needed in order that the VFT divergence at small x 
transforms smoothly to the Arrhenius law at large x- 
Previous calculations have used 6 = 3. 

The equation of motion for x itself is a statement of 
the first law of thermodynamics; it describes entropy flow 
through the slow, configurational degress of freedom into 
the fast thermal motions of the grains. Near steady state, 
it has the form 



X oc e 



-i/x 



r(x(<?) - x) + K P 



X 



(5.4) 



The first term in the square brackets on the right-hand 
side is the rate at which x is driven towards x(q) by 
the mechanical noise strength T. The second term, pro- 
portional to p, is the rate at which thermal fluctuations 
drive x toward the scaled ambient temperature 9. k is a 
dimensionless parameter of the order of unity. The com- 
petition between the two terms in Eq. (|5.4|) determines 
the value of x; it is close to xil) f° r large T, and close to 
9 when the system is driven slowly so that T is small. 

There is one complication that must be dealt with at 
this point. Equation (|5.4[) . as written, implies that the 
steady-state x must lie in the interval between xil) and 
9. If we assume that x(q) ~ Xo is a constant for small 
enough q, then a system initially prepared with 9 > xo 
would be "cooled" to x < $ when driven at a small strain 
rate. This behavior seems implausible; so far as we know, 
it is not seen in simulations, e.g. [14| . in which x 1S 
measured directly. In [l9T | , this problem was corrected by 
setting xo = when 9 exceeds xo> and by rescaling xi 
and xa accordingly. Specifically, 



Xo 



Xo, if < xo; 

e, n9>xo, 



(5.5) 



and 



Xi = 



XA 



Xi, if ( 
OXi/Xo, if' 



< xo; 
> Xo- 



XA, if < Xo; 

Q~Xa/xo, if > Xo- 



(5.6) 



This guarantees that x(q) > at all times. 



In summary, we use Eq. (|5.4p in the form 
T X (q) + K P 9 



X : 



r + np 



(5.7) 



along with Eq. (|5.2I) . to determine one relation between 
X, q and s. We then use Eq. (|4.6I) to compute both q and 
X as functions of s. 



VI. ANALYSIS OF THE HSL SIMULATIONS 

In their molecular dynamics simulations, HSL studied 
mixtures of 4096 hard spheres, half each of size a and 
1.4 a, in contact with a thermal reservoir at fixed tem- 
perature 9 = 1, undergoing steady-state, simple shear 
driven by a shear stress s and a pressure p. They re- 
ported the strain rates as functions of s/p for a wide 
range of dimensionless temperatures 



1H 



pa 



(6.1) 



and packing fractions <f>. In a separate paper [l4j ]. Haxton 
reported measurements of the compressibility tempera- 
ture 6* C omp of this system using the relation 



{N/V) K fl^p = S{k -> 0) 



(6.2) 



where K is the isothermal compressibility and S(k) is 
the structure factor as a function of the wavenumber k. 
He argued that 9 comp is the thermodynamic temperature 
of the configurational degrees of freedom for this system, 
and therefore ought to be the configurational effective 
temperature in nonequilibrium situations where 9 comp ^ 
9. Accordingly, in analogy to Eq. (|5.1|) . we assume that 



X 



pv z 



(6.3) 



A selection of five of the thirteen HSL data sets is 
shown in Figures 1-3, along with our corresponding 
theoretical results. Figures 1 and 2 are log-log plots, re- 
spectively, of q and the dimensionless, viscosity-like ratio 
s/pq as functions of s/p, for the values of 9h indicated 
in the figure caption. In Fig. 3 we show graphs of 



_ "comp _ Vz 

Xa = o r X 
pa 3 a A 



(6.4) 



as functions of log 10 (g), again for five selected values of 
Oh- 

The following analysis is closely parallel to that which 
has been used for interpreting molecular dynamics simu- 
lations of shear flow in a metallic glass [22[ ■ In [22[ , the 
system undergoes an ordinary thermal glass transition, 
but there were no direct measurements of the thermody- 
namic effective temperature with which to compare the 
theory. As will be seen here, the measurements of Xh 
shown in Fig. 3 provide crucially important information. 
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\og l0 (s/p) 



FIG. 1: (Color online) Log-log plots of simulated (data points) 
and theoretical (solid curves) values of the dimensionless 
strain rates q = ry as functions of s/p, for dimensionless tem- 
peratures 6h = 0.03 (blue circles), 0.05 (green squares), 0.08 
(orange diamonds), 0.10 (magenta triangles) and 0.25 (red 
down triangles), shown in that order from bottom to top. 



We find that we can fit the HSL data with a highly 
simplified version of the STZ theory outlined above in 
Sees. Ill - V. Our first major assumption is that the 
quantities p, sq, 7£(0) = i?o and q$ are functions only of 
the scaled temperature 9 defined in Eq. (|5.ip . 

Second, we assume that the symmetric function C{s) 
defined in Eq. (|3. 12|) is simply a constant, i.e. C(s) k, 
C(0) = Rq{9). We further assume that the inertial rate 
t -1 is a sufficiently accurate estimate of the attempt fre- 
quency for high-temperature, activated processes that we 
can let Ro(9) — > 1 and p(9) — > 1 in the limit of large 6. 
We do not attempt to predict the temperature depen- 
dence of either of these fundamentally glassy rate factors 
at low temperatures; rather, they are regarded here as 
measurable quantities. 

Our third physically significant assumption is that the 
ratio a 3 /vz is a constant, independent of 0. In other 
words, we assume that the activation volume relevant to 
STZ creation is some fixed multiple of the volume of a 
grain, independent of whether the system is above or be- 
low a jamming transition, independent of whether the 
transition rate Rq{0) is fast or slow, and even indepen- 
dent of whether the STZ's are compact or extended ob- 
jects under such different circumstances. Our idea is that 
a 3 is a measure of the free volume that must be created in 
order for grains to exchange places with each other, and 
that it does not matter whether that volume is created 
quickly in a small, loose cluster of grains, or whether it 
is created more slowly via rare fluctuations in a larger 
region. For similar reasons, we assume that 2vq, the vol- 
ume of the plastic core of an STZ, is also a fixed multiple 
of a 3 . 

Finally, the dimensionless numbers eo and k, which 
first appeared in Eqs. (|3.11[) and (|5.4I) . should both be of 
the order of unity. In the absence of better information, 
we set cq — n = 1. 




- n i 

-3.0 -2.5 -2.0 -1.5 -1.0 -0.5 
logi Cs/p) 

FIG. 2: (Color online) Log-log plots of simulated (data points) 
and theoretical (solid curves) values of the dimensionless vis- 
cosities s/p q as functions of s/p, for dimensionless temper- 
atures 8h = 0.03 (blue circles), 0.05 (green squares), 0.08 
(orange diamonds), 0.10 (magenta triangles) and 0.25 (red 
down triangles), shown in that order from top to bottom. 

The ratio vz/a 3 is especially important for our pur- 
poses, because it is the factor needed to convert from our 
temperature variables 9 and x to the HSL variables 9h 
and xh, as shown in Eqs. (|6.ip and (|6.4|) . We can deter- 
mine this ratio, as well as the value of Vq/vz, by looking 
at the HSL data in the limit of high temperatures and 
small stresses, where — as seen in Figs. 1 and 2 — the 
system exhibits simple linear viscosity. In this limit, we 
know from Eq. (|5.7|) that x ~ 9. We also know by taking 
the small-stress limit of Eqs. (|4.3p and (|4.6p . and setting 
p = Rq = 1, that 



Evaluating the left-hand side of this equation directly 
from the HSL data at temperatures 9h = 0.25 and 0.30, 
we find that vz/a 3 — 9h/9 = 0.53 and vo/a 3 = 1.65. 
Thus, the activation volume vz is somewhat smaller than 
the average grain size, and the plastic core of an STZ 
is somewhat bigger. Both are plausible results; but we 
emphasize that both depend on our assumption about 
the high-temperature limits of Ro and p. 

The next step in our analysis is to look at low temper- 
atures, in the jammed region, where p = and the data 
seem to indicate well defined yield stresses. In this limit, 
X = x(q), and Eqs. (|4.3I) and (|4.6I) become 

g^^Je- 1 ^) |l-^j. (6.6) 

This equation allows us to evaluate the function sq(6) 
directly from the observed yield stresses. More inter- 
estingly, Eq. (|5.2p tells us that x 1S a universal func- 
tion of the ratio q/q (9) so long as the parameters A, 
Xo, etc. appearing there and in Eq. (|5 . 3[) are indepen- 
dent of 9, which we assume to be true. Therefore, in 
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FIG. 3: (Color online) Log-linear plots of simulated (data 
points) and theoretical (solid curves) values of the effective 
temperature \h as functions of the dimensionless strain rate 
q, for dimensionless temperatures 9h = 0.03 (blue circles), 
0.05 (green squares), 0.08 (orange diamonds), 0.10 (magenta 
triangles) and 0.25 (red down triangles). 



this low-temperature regime, graphs of x(q) should col- 
lapse to a single curve when shifted along the log 10 (g) 
axis by amounts log 10 (<Zo)- The lower graphs shown in 
Fig. 3 do behave in roughly this way. Thus, by fit- 
ting these curves and using Eq. (|6.6I) to fit the low- 
temperature graphs in Figs. 1 and 2, we can determine 
the temperature- independent quantities A, xo, Xii an d 
Xa, defined in Eqs. (|5.5I) and (|5.6[) ; and we can deter- 
mine the low-temperature values of the functions Rq(0) 
and qo(9). Specifically, we find that A = 1.09, Xo — 0.12, 
Xi = 0.038, and xa = 0.47. 

Having evaluated the temperature-independent pa- 
rameters, it is a straghtforward exercise to complete this 
data analysis using the full STZ equations. We need 
to fit the data in the intermediate regions between fully 
jammed (p = 0) and fully unjammed (p — 1) situations. 
We also need to look in more detail at situations in which 
the driving forces s /p and the resulting strain rates q be- 
come large. In the latter cases, we know that the STZ 
theory must fail if the density of STZ's becomes too large 
and the system begins to behave more like a liquid than 
a solid. The HSL data shown in Fig. 1 do show signs of 
softening at the highest strain rates, and we make no at- 
tempt to fit such behavior. (In [22jj, it was found that this 
behavior at large stresses and low temperatures could be 
accounted for in part by addding a Bagnold-like stress 
dependence to the STZ transition rate.) On the other 
hand, the high-temperature data indicate nonlinear re- 
sponse at the larger values of q. This nonlincarity is sen- 
sitive to sq(9); thus we are able to evaluate that function 
outside the jammed region. 

The function p(6), shown in Fig. 4 as a function of 0h, 
indicates that the hard-sphere system undergoes a jam- 
ming transition in the region 0.05 < 6h < 0.10. Below 
those values of 9h , the rate factor p becomes unmeasure- 
ably small; but we have no way to determine whether 
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FIG. 4: (Color online) Thermal noise strength p as a function 
of the dimensionless temperature Oh- 



or not it actually vanishes. The internal transition rate 
Rq(0), shown as a function of 9h in Fig. 5, undergoes 
a qualitatively similar transition at smaller values of 6h, 
but remains nonzero down to the lowest temperatures for 
which data is available. Consistent with our discussion in 
Sec. V, the independently determined low-temperature 
values of qo(9) are roughly proportional to 4Rq(9), as 
shown in Fig. 5. This factor 4 might be the same fac- 
tor that appears in Eq. (16.61) ; but we emphasize that 
Eq. (|5.2I) . in which q Q was defined operationally, is no 
more than a phenomenological interpolation formula that 
should not be interpreted too literally. At higher temper- 
atures, where x ~ 9, qo(9) plays no role in the theory and 
therefore, in principle, cannot be evaluated. The stress 
scale so shown in Fig. 6 decreases as a function of 9r 
within the jammed region, and then becomes a constant 
at higher temperatures. This is qualitatively the same 
behavior that was observed in [19( and [22j. 

As shown in Figs. 1 and 2, the resulting fits to the 
data for strain rates as functions of stress are quite good 
everywhere except at the largest values of q where, as 
noted earlier, we expect the STZ theory to fail. The 
graphs of xh versus q in Fig. 3 are harder to decipher, 
in part because the data are noisier and, in part, be- 
cause the small-g behaviors happen to lie on top of one 
another. The two curves that are in the jammed region 
where p — 0, with 9h = 0.03 and 0.05, are slowly (log- 
arithmically) approaching xh — 0.53 xo — 0.064; but 
they are still above that value at the smallest q'a. The 
next two curves, for 9h = 0.08 and 0.10, already have 
reached their asymptotic values at xh — 9 Hi because 
r <C p in Eq. (|5.7p . As a result, the four lower curves 
in the Figure appear to be indistinguishable but, in fact, 
are consistent with the theory and reveal the importance 
of Haxton's effective-temperature data for understanding 
the competing time scales in this system. The single, dis- 
tinguishably high-temperature curve in Fig. 3 is nearly 
flat at xh = 9h = 0.25. The fact that this curve is 
crossed by the lowest-temperature curve with 9h — 0.03 
illustrates the sharp difference in the dynamics between 
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FIG. 5: (Color online) Dimensionless STZ transition rate 
Ro (blue circles) as a function of the dimensionless tempera- 
ture 9h- The red squares indicate independently determined 
values of go/4 at the four lowest values of 9h- Our rough 
estimates indicate that qo levels off at about 4.5 Ro above 
9h ~ 0.06, but is intrinsically not measureable at appreciably 
higher temperatures. 



the jammed and unjammed regions. 



VII. CONCLUDING REMARKS 

The thermalized hard-sphere model studied numeri- 
cally by HSL is an ideal laboratory for studying glass-like 
behavior near a jamming transition. It is sharply defined, 
with no ambiguities about the role of interaction ener- 
gies or about the origin of the underlying time scale. As 
shown in Sec. II, it lends itself easily to a nonequilibrium 
thermodynamic analysis that systematically constrains 
the STZ equations of motion and relates the Edwards 
compactivity X to the thermodynamic effective temper- 
ature x- 

The HSL data is especially revealing because of Hax- 
ton's direct measurements of the effective temperature 
[Tij , without which we could not have sorted out the tem- 
perature dependences of the different rate factors that 
appear in the STZ theory. For example, there was no 
such information about \ m the metallic glass simula- 
tions used by Egami and JSL [22|], and therefore no rea- 
son to expect that the strain-rate parameter qo defined 
in Sec. V was a strongly varying function of temperature. 
As a result, it was concluded there that the discrepan- 
cies between the theoretical and simulated viscosities at 
small strain rates and temperatures near the glass tran- 
sition were due to a failure of the single-species STZ the- 
ory. Here we see no such strong discrepancy, presumably 
because we are taking more accurate account of the 6 
dependence of the STZ transition rates in both Rq(9) 
and qo(9). This speculation, however, leads to yet more 
interesting questions. 

The most important of these questions, in our opin- 
ion, pertain to the relations between the present results 
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FIG. 6: (Color online) Stress scale so in units of 0/a A as a 
function of the dimensionless temperature 9b ■ 



and the postulate in [12| that the STZ's are the dynamic 
heterogeneities frequently invoked as explanations for a 
wide range of glassy phenomena. In [l(J EH, Bouch- 
binder and JSL showed that both basic statistical prin- 
ciples and the results of frequency-dependent viscoelas- 
tic measurements implied the existence of a broad spec- 
trum of STZ transition rates in equilibrated systems near 
their glass temperatures. In [12| . it was argued that 
the "slow" STZ's near the bottom of this multi-species 
STZ spectrum are responsible for anomalously large vis- 
cosities consistent with observed Stokes-Einstein viola- 
tions; and that the extended spectrum naturally pro- 
duces stretched-exponential relaxation in a variety of ex- 
perimental situations. 

Here, we have used only a single-species STZ theory; 
but we can interpret our values of Ro{9) as weighted av- 
erages over the full spectrum in the multi-species theory. 
A difficulty with this interpretation is that, according to 
the earlier statistical analysis, the spectrum should shift 
to higher transition rates as the system becomes noisier 
at larger driving stresses; but, so far, we have assumed 
that Ro(6) is independent of the stress. Nevertheless, 
the present results suggest that the statistically effective 
STZ's become markedly slower as 6 decreases through 
the jamming transition. Recently, Durian and coworkers 
[23l l24j. in experiments using air-fluidized grains, have 
found that their dynamical heterogeneities become slower 
as the jamming transition is approached. They also have 
found that these slower heterogeneities involve motions 
of larger numbers of grains, a size effect that seems es- 
pecially intriguing. It was argued in [l2j that the length 
of the diffusion jumps associated with slow STZ's must 
grow with decreasing temperature in order to account 
for the apparently normal, Fickian diffusion observed in 
those conditions. It will be important to find out whether 
these are related phenomena. 
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